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Abstract In this work we analyze in detail the behaviour and properties of the kinks 
found in an one dimensional model for the close packed rows of potassium ions in 
mica muscovite. The model includes realistic potentials obtained from the physics 
of the problem, ion bombardment experiments and molecular dynamics fitted to 
experiments. These kinks are supersonic and have an unique velocity and energy. 
They are ultradiscrete involving the translation of an interstitial ion, which is the 
reason they are called crowdions. Their energy is below the most probable source of 
energy, the decay of the isotope and above the energy needed to eject an atom 
from the mineral, a phenomenon that has been observed experimentally 
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1 Introduction 

Some materials are able to record the passage of charged particles and are used as 
radiation detectors ||2l|9l and there are minerals that show nuclear tracks that were 
produced at some stage during their formation ||22|. The mineral mica muscovite 
has been shown to have recorded the tracks of muons, positrons and other swift 
particles with positive charge while being deep underground Il23l425l . The most re¬ 
cent reviews are Refs. ll^lZTl . The tracks are recorded within the cation layer 
of potassium ions which form a two-dimensional hexagonal lattice. There are also 
many tracks along the close packed directions of this lattice that cannot be produced 
by charged particles and are attributed to some vibrational entities called quodons 
because of their quasi one-dimensional structure Il2^29ll3ll[^ . Their existence has 
also been shown directly with an experiment in which the energy of alpha particles 
incident on one side of a mica specimen was able to eject atoms at the opposite 
border along the cation lattice directions ll30ll . 

Recently, a model with realistic potentials for the dynamics of potassium ions 
within the cation layer of mica muscovite has been introduced |[T}{3l. The au¬ 
thors have considered the available potentials for the interaction between atoms and 
ions. For the interaction between potassium ions K+ the electrical potential was 
not enough because the passage of the kink brings about very short distances, for 
which the ions can no longer be described as point charges. Therefore, the Ziegler- 
Biersack-Littmark (ZBL) potential was used ll^ . This potential models the electri¬ 
cal repulsion by the nuclei partially screened by the electron cloud. ZBL potentials 



Fig. 1 Energies of several kinks with respect to time. When more energy than the crowdion’s one is 
delivered and therefore a faster kink is produced, a radiation process takes place until the supersonic 
crowdion is formed. Thereafter, the crowdion is extremely stable. If the initial energy is smaller 
than the crowdion’s one the kink dissipates into phonons. The scaled units are approximately 3 eV 
for energy and 0.2 ps for time. The final velocity and energy are approached asymptotically, being 
Vc = 2.7387 (7.2 km/s) and = 9.5 (26.2 eV) in scaled and physical units 
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have been widely tested and compared to data obtained in ion bombardment ex¬ 
periments, being therefore the more realistic ones while using classical mechanics. 
Quantum calculations could certainly provide more accuracy but at the cost of much 
more complex analytical and numerical calculations. The interaction of the potas¬ 
sium ions with the lattice was described with empirical potentials used in molecular 
dynamics and fitted with thermodynamic properties, neutron 0 and infrared spec¬ 
troscopy 0 and also validated for other silicates ifTTll . 

Arguably, the most important result in the full system with substrate was that a 
supersonic kink was formed with specific energy and velocity 0 . As it involves the 
movement of an interstitial atom through the lattice, it will be called a (supersonic) 
crowdion in this work as described in Ref. 03. The term will be reserved for this 
specific supersonic kink with stable and unique velocity and not for other kinks. If 
the lattice was given more energy, nonlinear waves and later phonons were emitted 
until the specific velocity and energy was reached. This characteristic of supersonic 
kinks associated with specific values of the velocity have also been described in 
Refs. 032II37L where they use the terms topological soliton and lattice soliton. 

The supersonic crowdion found in the mica model is extremely discrete as ba¬ 
sically only two ions are moving at the same time, which will be referred to as 
the magic mode with sinusoidal waveform and corresponds to a phase delay very 
close io q = 27r/3 II14II15I as explained below. In the magic mode, which was in¬ 
troduced in the Fermi-Pasta-Ulam lattice to describe both steady-state or slowly- 
moving breathers and supersonic kinks m, only two particles are mostly involved 
in the motion at the same time. The mode with mode q — 7t is the limit of discrete¬ 
ness as only one particle is moving at the same time, and the kink is equivalent to just 
one particle hitting the following one with a hard-sphere interaction. We have also 
called these kinks ultradiscrete kinks (UDK). They are also known as kinks with 


Fig. 2 The structure of mica 
muscovite where the potas¬ 
sium layer can be observed. 
This point of view has been 
chosen to emphasize the K+ 
rows represented by yellow 
balls 
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atomic scale localization and have been described theoretically IfTOll and observed 
experimentally in a chain of repelling magnets ES). The energy dissipated by the 
crowdion and its subsequent stability can be seen in Fig. [T] Supersonic kinks with 
a discrete set of velocities for which there is no radiation have been described in 
previous publications 013lll8ll^lJ7l . They appear in systems with substrate poten¬ 
tial and nonlinear coupling and can be described as multiple solitons. In our system 
due to the extreme discreteness of the kinks there is only a non-radiating velocity 
corresponding to a double soliton as will be explained in Sect.|4] See also Ref. 13. 
Figure 13 represents the coordinates of the potassium ions obtained in a numerical 
simulation. 

The energy of the crowdion is approximately Ej^ = 26.2 eV, which is an interest¬ 
ing result because there are sources of energy in the lattice which can provide it as 
it will be explained with more detail in Sect. |7] The most abundant of the unstable 
potassium isotopes is which can decay by different beta processes providing 
recoil energies up to 50 eV. The crowdion energy is also smaller than the second 
ionization energy of potassium, that is, 31,6eV El, which thus prevents the pos¬ 
sibility of inelastic collisions where the kinetic energy would be lost stopping the 
propagation of the kink. It is also larger that the 4-8 eV needed to eject an atom iTf^ . 
an effect that has been found in an experiment where the transmission of localized 
energy along lattice directions with the subsequent ejection of an atom at the edge 
of the boundary has been observed ll^ . 

Another point of interest of the crowdion is that it consists essentially of a 
charged interstitial K^, i.e. an excess of an unit of elemental charge, travelling at 



t 

Fig. 3 Coordinates of the supersonic crowdion or ultradiscrete kink from numerical simulations. 
It can be observed that only two particles are moving at the same time. Lattice units a — 5. 19A for 
coordinates and scaled units (0.2 ps) for time. Also the double kink structure can be seen as will be 
explained later in the text 
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twice the speed of sound. Therefore, it is very likely to be recorded, as positively 
charged particles leave tracks in mica muscovite. 

Are the quodons observed in mica muscovite the crowdions described in this 
work? It is not clear, but there are several points in their favour; a) They have an 
energy that can be produced by the recoil of b) They have enough energy 
to expel an atom at the surface; c) They have stability and seem to travel forever; 
d) They survive to room and higher temperatures; e) They transport positive charge 
that would leave a track in mica muscovite. Against them is that their existence and 
stability has not been verified in two or three dimensions. But even if their energy 
spreads they are likely to leave some of the other dark marks in mica. 

The sketch of this work is as following; Section |2] describes the system and po¬ 
tentials. In Sect.[3]the magic mode is described with detail and the quantities in the 
fundamental ansatz are redefined in a new meaningful way. Section|4]describes the 
properties of the kinks when the substrate potential is introduced and the supersonic 
crowdion appears, while Sect.|5]describes the properties of phonons in a system with 
a substrate and applies them to analyze the crowdion’s phonon tail. Some interest¬ 
ing results of the outcome of numerical simulations when excess energy is delivered 
and when the system is previously thermalized are presented in Sect. |6] The recoil 
energies in the different decay modes of and their consequences for the forma¬ 
tion of kinks or other lattice excitations are described in Sect.|7] The work ends with 
a summary. 


2 Description of the system 

Mica muscovite is a layered silicate where a layer of potassium ions is sandwiched 
between layers of a complex silicate structure of tetrahedra and octahedra. This 
cation layer has a hexagonal structure where rows of potassium ions can be iden¬ 
tified, as seen in Fig. |2] As explained with more detail in Refs. flJQ we consider 
an ID model for a row of K+ ions. The distance between ions is a = 5.19A which 
in scaled units will be taken as the unit of distance. The interaction between ions is 
described by two terms, the first one is the electrostatic Coulomb repulsion 

2 2 

Uc = Ke--Ke-, ( 1 ) 

r a 

where Kg is the Coulomb constant, e the elementary unit of charge and r = 
d„ = Xfi — Xn-\ is the interatomic distance. The reference level of energy is taken 
as the electrostatic energy at the equilibrium distance a. This value of energy 
Kg^ja = 2.11 A6eW is also taken as the unit of energy in scaled units, and it is 
useful to remember that it is approximately ue ^ 3eV. The other natural units 
are the potassium mass iuk = 39.1 amu and therefore the derived unit of time 
T = \/mKa^/Kgt^ — 0.1984ps ~ 0.2ps. 

This system supports propagating kinks of almost any velocity and energy IMl 
but with very small inter-particle distances for which the ions cannot be described as 
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point particles. The second term for short-range repulsion is the Ziegler-Biersack- 
Littmark or ZBL potential, which corresponds to the electrostatic interaction be¬ 
tween nuclei partially shielded by the electron cloud which is described by an uni¬ 
versal function that has been tested with experiments of ion bombardment ll^ . The 
ZBL potential usually involves four terms which are effective at different ranges of 
energies. For the potassium atoms at energies up to 200keV it is enough to consider 
a single term given by 

t/zBL(0 =-exp(--^), (2) 

r p 

with a = 2650.6 eVA and p = 0.29529 A which correspond to a = 184.1 and p = 
0.0569 in scaled units, respectively. No attractive Van der Waals term is considered 
as it would be much weaker than the repulsive term. The system with Coulomb and 
ZBL potential also support propagating kinks with many energies but with realistic 
distances between particles 0. 

The interaction with the atoms in the lattice above and below the potassium layer 
is obtained from an unrelaxed lattice using empirical potentials used in molecular 
dynamics and fitted with thermodynamic and spectroscopic properties Il5]fm which 
are also valid for other silicates. The resulting periodic potential can be written as a 
Fourier series for which it is enough to retain the first five terms 0 

Us{x) = 52 t4cos(27rn-). (3) 

«=0 



Fig. 4 Interaction potentials U(r) in scaled units. Coulomb (—); ZBL (-); Coulomb+ZBL (thick 

-); substrate potential (• • •) and the sum of the Coulomb, ZBL and substrate potentials (-). The 

scaled units are 2.77 eV and the lattice unit a = 5.19A for U and r, respectively 
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The Fourier coefficients are given by 

U„ = [6.7902,-9.2920,3.0512,-0.6387,0.0891]eV = 

= [2.4473,-3.3490,1.0997,-0.2302,0.0321], (4) 

with the latter values given in scaled units. As will be shown later, the linear spatial 
frequency for the long wavelength limit becomes 119 cm^', that is quite close to the 
experimental one of 110 cm^' obtained with infrared spectroscopy m. A compari¬ 
son between the three potentials can be seen in Fig. |4| 


3 The magic mode revisited 

In this section we describe the fundamental ansatz and the variable involved. We will 
define the variables in a proper way, as they are not the same as in plane waves in 
spite of their analytic similarity. We will use scaled units for which the equilibrium 
interatomic distance is the unity as described above except where stated otherwise. 


3.1 Basic variables 

Some variables used throughout the study are introduced here, together with their 

definitions: 

Position .r„: It describes the position of the particle labelled n. At equilibrium 
Xn = na, although the origin of n is arbitrary. 

Displacement u„: It is the separation of the particle n from the equilibrium posi¬ 
tion, that is u„ = x„ — na. 

Interatomic distance dn'. It is the distance between two particles or ions. At equi¬ 
librium it is equal to the lattice unit a, which in lattice units is the unity, but 
it will be written often explicitly for clarity. It is related with the positions and 
displacements as d„ = Xn— x„_ i = u„ — m„_ i -f a. 

Strain The increase of d„ with respect to the equilibrium distance, i.e. v„ = 
dn — a. It is always negative for the kinks described in this work. It is related with 
the displacements as: y„ = — m„_i. 

Compression c„: The decrease of dn with respect to the equilibrium distance, i.e. 
Cn— a — d„ = —Vn- It is always positive for the kinks described in this work. It is 
related to the displacements as: c„ = 
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3.2 Fundamental ansatz 

As demonstrated in Refs. OlIISl for a large set of kink solutions of Fermi-Pasta- 
Ulam systems, the strain v„ = i can be approximately described by the 

fundamental ansatz with sinusoidal waveform: 

Vn =— — {1+cos{q{na — Vt)) with — n < q{na — Vt) < n, (5) 

where q = Injlia or q = 2n/l> in scaled units with a = 1 that we will usually use. 
The value of v„ is always negative representing a compression of the bond. This 
ansatz describes a moving profile with velocity V that it is better visualized in the 
alternative form v„ = —Acos^{q/2{n — Vt)). At any given time its value is zero 
except for a length A =2nlq representing the number of consecutive bonds com¬ 
pressed. For a given bond n the value of v„ is zero except for an interval of time 
T = 27il{qV) representing the time during which the bond is compressed. Note that 
A is not a wavelength as there is no periodic wave and T is not a period as there is 
no periodicity in time. 

For convenience we will often use the equivalent expression for the compressions 

Cfi = Vn'‘ 


c„ = —(1+cos(a)f — ^n)) with —K<COt — qn<K, (6) 

where co = qV is the rate of variation of the phase ^{nf) — cot — qn, i.e., co = 
d^{n,t)/dt but it is not the frequency as there is no periodicity. This equation will 
be used in the next subsection as it is easier to interpret because c„ is always positive, 
the phase increases in time and the bonds compressed later have smaller phase. 

From the fundamental ansatz the displacement can be constructed and it may be 
instructive to compare them with other solutions. They can be seen in Fig. |5]for 
the magic mode q = 2n 13 compared with the first solutions for supersonic crow- 
dions ns. The compressions c„ = have a solitonic form and in the same 

figure they are compared with the discretization of the solutions for the KdV equa¬ 
tion, which describes waves in a canal m, one of the first examples of solitons. 

As we have seen these equations are not as simple as they seem, due to the com¬ 
pactness condition for being nonzero. They look like harmonic waves, but they are 
not. The quantities in the equation have to be redefined but they keep the usual re¬ 
lationships for harmonic waves. In the following we propose operational definitions 
that are convenient but are only approximately correct, which is also natural at the 
fundamental ansatz is not exact either. 

Velocity V: The average velocity of the kink. This is the magnitude best defined 
in numerical simulations and experiments. 

Phase (j) {n,t): Trivially, the phase of the bond nis^{nf) = (Ot — qn. It determines 
when a bond is compressed —n < ^{nf) < n and its state of compression. For 
example, (pinf) = 0 is the phase of the state of maximum compression of the 
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Fig. 5 (Left ) Profiles of the displacements u„ for the sinusoidal magic mode (o) and the ones 
given in the original supersonic crowdions paper by Kosevich and Kovalev (1973) (H. For a 
quartic interatomic potential (A): u„ — (2/;r)arctan[exp(—— W))] and for a cubic one (V): 
«;i = [1 + exp(2g(« — . ( Right) Comparison of the compressions c„(f) = «„-i — u„ for the 

magic mode (-) with the soliton for the continuous KdV equation fTH : c„ =Asech^[g{«-yf)] 
(- -). The functions have been rewritten so that the parameters have the same meaning. The magic 
mode is between the two K&K solutions and it is wider than the KdV one 


bond Cn =A, (j){n,t) = —;r means the beginning of the compression process and 
^{n,t) = n is the end. It is not periodic as a bond is just compressed once, if for 
example, (/)„ = 2;r c„ = 0 and not A. 

Active: This term will change depending on the variable we refer to. For the phase 
it corresponds to 0(n,f) S \—n,n). 

Phase rate (0: It is the rate of variation of the phase with time or (a — d(j){n,t)/dt = 
qV. It is not the angular frequency as the ansatz is not a periodic function. 

Compression time T: It is the interval of time for which a bond is compressed or 
activated, T — Injo). The interval of activity starts with zero compression c„=Q 
and finishes with the same value. In the meantime it achieves c„ = A, its maxi¬ 
mum value. It also starts with (j) = —7t and finishes with (j) = n. As the numerical 
solutions become separate from the fundamental ansatz the operational definition 
of T is the value that brings about a better fit of v„ with the fundamental ansatz. 

Phase delay q: It is the phase difference between two active (compressed) bonds 
n and « — 1, that is, q — (j){n,t) — ^{n — l,f). Alternatively, it can be defined as 
q — 2n{dt/T) — codt, where 5t is the time delay between two consecutive active 
bonds. 

Kink length A; It is the spatial extension of the kink, very much related with the 
number of active bonds at a given time X/a or simply X in scaled units. It is 
given by A =2%!q and it is also the distance travelled by the kink during a time 
interval T, i.e., A = VT. The usual relationships also hold, that is, V = (o/q and 
A = 2%! q. 

Amplitude A: It is the maximum value of the compression c„. 
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Fig. 6 Fourier spectrum of 
the kinetic energy of the su¬ 
personic crowdion obtained 
from numerical simulations. 
It is measured in a frame 
that moves with the crowdion 
in the lattice. We use arbi¬ 
trary units for the intensity 
and scaled units (5 THz) for 
the frequency. The value of 
the first harmonic is exactly 
the characteristic linear fre¬ 
quency V = Vela = 2.7387 
and circular frequency 
S) = 2;rv ~ 17.2, which 
corresponds to v ~ 13.4THz 
in physical units. 


Minimum distance R: It is the minimum value of the interparticle distance dn, that 
IS, R — a — A or R = 1 — A in scaled units. 

Characteristic frequency V: This is the inverse of the time 5t that the kink needs 
to travel a distance of a lattice site, i.e. v = 1 /dt — V/aor simply v = V in scaled 
units. Note that v = {X/a){\/T) (and not 1 /T). As the kink is not periodic it is 
the physical frequency at which the compression, the kinetic or potential energy 
or other magnitudes change while the kink travels in a lattice with period a. 
An example can be seen in Fig. |6] Their values for the crowdion are therefore 
V = 2.7387 and d) = 2;rv ~ 17.2, corresponding to v ~ 13.4 THz. 

The equations for the displacements m„ and its derivatives will be seen in the 
following subsection. 


3.3 Phasors for the magic mode 

The easiest way to visualize the relative phases and distances of the variables is to 
consider the rotating complex vectors or phasors 

b„ = with ^{n,t) = (Ot — qn and c„ = ^-f Re(b„), (7) 

There is an important difference with the usual concept of phasors and it is that 
the circle is not periodic. The only phase interval where the phasors exists is — tt < 
(j){n,t) < n. If (j){n,t) < —7t the phasorb« has not yet come into existence and when 
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^{n,t) > n, b„ has disappeared. Therefore, for q = 27r/3 at a given time there are 
three phasors in the unit circle as shown in Fig. |7] The three phasors have their 
origin at (A/2,0) and rotate anti-clockwise with angular speed (O while the time t 
increases, let us denote them b„_i, b„, b„+i. In the following n has to be understood 
as the index of the inner bond of the three compressed ones or the index of the 
intermediate phasor, that is —;r/3 < ()>(«,f) < 7r/3. If we denote as = n/V the 
time for which = 0, then ~T I6<t — t„ < T/6. This is not a restriction as 

there is always a bond central to the three compressed ones. 

The phasor b„+i is behind b„ ny an angle q and so on for a kink travelling to 
increasing n number. Note that b„_i + b„ + b„+i = 0. 

Therefore, the particles first reached by the kink have larger phase The angle 
^ = TT is the angle for change of number, that is, when b„_i reaches that position 
it disappears from the circle and ceases to be active, indicating that the bond n—\ 
is no longer compressed. At the same time, a new phasor b „+2 appears at 0 = —7t, 
indicating that a new bond has started to be compressed or becomes active, after a 
time T it will in turn become inactive. As shown in Fig.|2l the horizontal distance to 
the vertical straight line through the origin is the compression c„ = A/2 + Re(b„). 

Let us now consider the displacements u„, using c„ = m„_i — m„ or m„_i = + 

Cn- The particles not yet reached by the kink have zero displacement and the first 
nonzero compression is c„+i. Therefore = Cn+i and —Un + c„= c„ +c„+i = 



Re 

Fig. 7 Visualization of the evolution of the compressions c„ = —v„ = m„_i —u„ for the magic 
mode q = 27t/3 (A = 2/3). Three phasors b„_i,b„,b„+i centered at (0,A/2) and rotating anti¬ 
clockwise are active (the bonds are compressed) at a given time t when —n < (j){n,t) < 7t. Their 
horizontal coordinates give the compression as can be seen with c„+ 1 . The maximum compression 
A is achieved for (j>{n,t) = 0. At ^(n,t) — it, h„ i will transform into b „+2 indicating that the bond 
n — I is no longer compressed while the bond n -|- 1 starts its compression cycle. The displacements 
are active while changing and only two are active at a given time u„ = c„+i and m„_i = c„ + 
Cn+i —3Aj2 — c„-i. For m > n, = 0 and for m < « — 1, m„, = 1. Also the nonzero velocities are 
u„ = —(Blm(b„+i) and = t()Im(b„-i). Magnitudes are in lattice units a = 5.19 A 
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A + Re(b„ + b„+i) = A + Re(—b„_i) = 3A/2 — c„_i as represented in Fig. [7] To 
summarize 


Un-\-l — 0 ( 8 ) 

A A A 

= Cn +1 = - +Re(b„+i) = - + -cos{cot-q) 

3A A 

Un-i = —-Cn-i =A-Re(b„_i) = A - -cos{cot + q). (9) 

These equations are valid for t — 0 chosen as the time for which the central bond 
n is most compressed Cn = A and remains central, —7zj3 < (j){n,t) < n/3 and 
—T l6<t<Tl6. The following displacement m „_2 = c„-i +c„ +c„+i — 3A/2 — \ 


-3T/6 

-2T/6 

-T/6 

0 

T/6 

2T/6 

3T/6 


n-3 n-2 n-1 n n+1 n+2 n+3 

Site 

Fig. 8 Magic mode q — iTz/'i for a kink. A sketch of the system is shown for a full time of 
compression T at time intervals 7/6. The white particle is labelled n, therefore its displacement 
is u„ and the bond at its left is also bond n with compression c„ = — v„ = m„_i — u„. The origin 
of time has been taken as the time of maximum compression of bond «, i.e., c„ = A and d„ — 
a — A. During the time interval in the graph the white particle n moves from site n to site n+1. 
At time t — —T /2 = —37/6, the bond n is uncompressed (c„ = 0, d„ = 1) and again becomes 
uncompressed att = 3T/6 = T /2. Note that during the first two 7/6 intervals, although the bond 
d„ is changing, there is no appreciable displacement u„. Note also, that the compressed structure 
at f = —37/6 between sites n — 2, n — 1 has moved at t = 37/6 to sites n + 1, n + 2, i.e., the kink 
has moved three sites or the length of the kink X =VcT, while the white particle n has moved a 
single site. Therefore, the average velocity of a particle in a time 7 is (V/,) = 1/7 = V'c/3. The 
average velocity of a particle for the following four 7/6 intervals, when it is actually moving, is 
(+„>'= l/(27/3)=T,/2 
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and equally = 1 for m <n —I, that is, for the particles that have been left dis¬ 
placed by a lattice unit after the passage of the kink. 

The particle velocities m,„ = dum/dt can also be calculated and visualized easily 
using b„, = icobm and therefore Re(bm) = —a)Im(bm) 

Un = Cn+\ = —a)Im(b„+i) = 2 “ ?) 

= —c„_i = ®Im(b„_i) = 01 —sin(a)f-l-g'). (10) 


For any other m, iim = 0. 

For other integer values of X =2nlq, there are X active phasors and for non 
integer values, the number of active phasors changes between the two integers below 
and above X. However, in this work we will concentrate on the magic mode q — 
27r/3 as it is very close to the crowdion found in the simulations. 

In this way it is easy to construct the evolution of the particles during the com¬ 
pression time T as can be seen in Fig.|8]for six times between —T /2 to T /2. In this 
time the crowdion advances a length A = 3, that is, three lattice units, but a single 
particle just travels a single lattice unit. Therefore the average velocity of a particle 
< Vp > is three times smaller than the crowdion velocity 14. It is worth mentioning 
that Fig. [8] also shows that only the two particles participating in the kink motion 
are mostly involved in the motion at the same time, as the fundamental ansatz with 
sinusoidal waveform, Eq.j^with q' = 27r/3 predicts 01411151 . 


4 Kinks with substrate potential: the crowdion 

The introduction of a substrate potential also modifies substantially the behaviour of 
the particles in the kink. The phase (j) {n,t) is still very useful for the interpretation of 
the movement of the particles. The crowdion, of ultradiscrete kink of fixed velocity 
and energy that appears in the simulations corresponds basically to the magic mode 
but with some differences. Considering the white ball in Fig.|8]and denoting it by n, 
it basically does not move from t G [—37’ jb^—T /6] as the Coulomb repulsion from 
particle n — 1 is weak. For times close to f = 0 when the strong ZBL potential acts, 
it receives most of its momentum which it will transfer in due course to the follow¬ 
ing particle n + \. However, in between, it will have to overcome the barrier of the 
potential, experiencing a deceleration and afterwards an acceleration while going 
downhill. Eventually the acceleration becomes negative as it experiences the ZBL 
repulsion from the particle n -f 1 ahead. The ascending and descending of the poten¬ 
tial barrier by the particle produces a remarkable change in the particle displacement 
M„ and strain v„ = u„ — m„+i as shown in Figs.|9] The kink has been converted into 
a double kink: the first kink corresponds to the translation of a particle from the 
well bottom to the top of the nearest potential barrier and the second kink to the 
subsequent displacement to the following well bottom. 
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t /T 

Fig. 9 Comparison of the ultradiscrete kink defined with the fundamental ansatz in Eq.|^with A = 
2/3 and q —271/3 with the ultradiscrete kink with fixed velocity obtained in the simulations dubbed 
crowdion in this work. Dashed lines: ansatz, continuous lines: crowdion. The displacements u„ 
correspond to the upper curves while the strains v„ — u„ — u„-i correspond to the lower curve. 
The kink transforms into a double kink because the displacement between two equilibrium sites 
is divided by the nonequilibrium position at the top of the potential well as can be seen in Fig. 
The magnitudes u„ and v„ are given in lattice units a = 5.19A. The compression time is given by 
T = 1.095 or 0.22 ps in physical units 


We would also like to mention in connection with Fig. |9l that the fundamental 
ansatz with sinusoidal waveform, Eq. (5) for q — In/h, and corresponding dashed 
lines in these hgures give much better agreement with the simulations of supersonic 
kink motion in the Fermi-Pasta-Ulam lattice without substrate US. The deviation 
from the ansatz prediction in Fig.|9]is caused only by the presence of the substrate 
because the ansatz was originally proposed for the translationally-invariant Fermi- 
Pasta-Ulam lattice iHiia. 

The separation from the ideal functions of the ansatz can also be seen in Fig. [TO] 
where the displacements are shown at a given time. It can be observed that the de¬ 
viation from the magic mode are important qualitatively but not so much quantita¬ 
tively. A more signihcant difference appears in the velocities which are represented 
in Ref. |0 but can also be seen easily in the slope of Fig.|9] According to Eq. (fTOl l 
the maximum particle velocity using the ansatz is (£)A/2 =1.91 or 5km/s, while for 
the observed one for the crowdion it is 2.9 in scaled units or 7.6km/s attained when 
the particle is going uphill or downhill. The minimum particle velocity is achieved 
at the top of the barrier. 
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Fig. 10 Three plots at different times (a) f ~ — r/6 (b) f ~ —0.5T /6 (c) f ~ 0. They show the profile 
of the displacements u„ (upper curves) and strains v„ — u„ — u„^i (lower curves) with respect to 
the particle and bond index n. The continuous lines represent the theoretical ansatz Eq. 0 ) and the 
circles represent the points corresponding to the numerical simulation of the crowdion. Time t — 0 
corresponds to the maximum compression of bond n. The variables u„ and v„ are given in lattice 
units a = 5.19 A. Every T/6 the theoretical and numerical solutions becomes almost identical as 
can also be seen in Fig. |9] Subfigure (b) represents the maximum separation from the theoretical 
curves 


5 Phonons and crowdions 

The introduction of the substrate potential brings about significant changes in the 
system, not only for the kinks but also for the phonon spectrum. We first review 
the properties of phonons in a system with substrate potential and then use them to 
analyze the phonon tail of the crowdion. 


5.1 Phonons in presence of a substrate potential 

The dynamical equations for small perturbations become 

Ufi — dip Ufi T C^. 1 2 Un) , (11) 

with Cs = \/2. The linearization of the coupling terms come only from the Coulomb 
one. The ZBL potential does not appear because it is zero for small oscillations. The 
substrate potential has been reduced to a harmonic one expanding the sinusoidal 
functions. The value of (Oq is obtained using the values of the Fourier coefficients of 
the substrate potential in Eq. (|4]i 

4 

®0 = - ^ {27tmfUm ■ 

m=l 


(12) 
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Fig. 11 (a) Dispersion relation, (b) phase velocity and (c) group velocity. The three plots are for 
longitudinal phonons in a potassium row for the system without substrate (dashed line) and with 
substrate (continuous line). Scaled units are approximately 5THz for frequency and 2.6km/s for 
velocities. White circles are measurements from different numerical simulations. The black circles 
are the theoretical values for the phonon tail obtained by making the phonon phase velocity equal 
to the crowdion velocity 


The resulting numerical value is Oo = 4.48 in scaled units, corresponding to 3.6 THz 
or 119cm^^ The coefficient Cs = v/2 or 3.7km/s in physical units is the speed of 
sound in the system without substrate. 

Substitution of m„ = exp(i(q'n — COt)) leads to 

-m2 = _®2_^c?(ei? + e“i^-2). (13) 

From this equation it is easy to obtain the phonon spectrum, the phonon velocities 
and the group phonon velocities. They are given by 

COq+ 4c?sin^(|) ; Vphase = — 

z q 

to ^ casing 

Y^®o+4c2sin2(|) 

The corresponding equations for the system without substrate are identical with 
©0 = 0. In this case c* is both the phase and group velocity in the long-wavelength 
limit (q ^0). 

For the system with substrate (Oq is the lowest phonon frequency, corresponding 
to the long wavelength limit (q —0). This can be seen in Fig.[TT]where the disper¬ 
sion relation, the phase and the group velocities are shown. Note the main changes 
produced by the introduction of the substrate potential: (a) the phonon spectrum 
becomes optical, i.e., bounded from below, (b) the phase velocity diverges when 
q ^ 0, and (c) the group velocity becomes zero both at q = 0 and q — n and has a 
maximum close to q — Kj2 but with a much lower velocity. 

The value of the wavevector q corresponding to the maximum group velocity 
can be calculated as it corresponds to dVg/dq — 0. Equivalently it corresponds to 
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the maximum of the function 

f(q) = ^ = __, 

cUfl+2c2(l — cos(^)) ’ 

where we have used that 2sin^(^/2) = 1 — cos(^). Then 

df{q) 2sin(^)cos(^)[(UQ +2c^ — 2c^cos(^)] — sin^(^)[2c^sin(^)] 
d^ (Oq+ 2cj— 2c]cos{q) 

Making the numerator equal to zero, we obtain; 

(oIq + 2c^) cos (q) — 2c^ cos^ (q) — sin^ (q) =0, 

which leads to a second order equation in cos(q) 

c^cos^(q) — (o)q +2c^)cos(q) + c^ = 0, 

with solution 

®o+2c?± J< + 4®2 c 2 
cos(^) =- ^2 -• (19) 

For the values in the present system, only the minus sign gives a real value of ^ = 
1.4870 rad corresponding to a wavelength A = 4.2253, and maximum group velocity 
yg,M = 0.4091. 


5.2 Crowdion phonon tail 

When the kink is produced, its amplitude diminishes towards the crowdion’s one in 
an asymptotic way. Therefore after some time, the nonlinear waves are no longer 
produced but there is always a linear vibration left behind although with decreas- 
ingly smaller amplitude. This is why the crowdion continues propagating. The tail 
is a plane wave and as such does not transport energy, but theoretically could be 
measured to detect crowdion properties. We will call it the phonon tail. Note that 
the velocity to describe these plane waves is the phase velocity and which in this 
case is unbounded. The crowdion is moving at speed Vc and leaves at each site some 
small perturbation exactly at the same estate at times separated by = 1 /Vc. In 
other words, the phase velocity of the phonon tail V is the same as the velocity of 
the crowdion Vc . 

Fphase = Vc = 2.7387 (7.2km/s) (Phonon tail). (20) 

The wave number of the tail can be obtained from the equation Vphase = Vc = 
w/q = [cOq +4sin^(^/2)]*/^/^, which can be solved numerically or graphically 


(15) 

(16) 

(17) 

(18) 
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Fig. 12 (Left) Plot of c„ + n where the double soliton structure, period and other crowdion param¬ 
eters can be appreciated. (Right) Phonon tail amplified 50 times. It is a perfect plane wave with 
parameters with the same velocity of the crowdion V = Vc and similar parameters although not 
identical T ^ A > Ac. Note that this parameters are not well defined for the crowdion 

and only approximate 


from Fig. fTTT b). The result q= 1.8290 = 0.5822;r and therefore co = qVc = 5.00, 
T = Injd) = 1.2544 and X =2nlq = 3.44. So the parameters are very close to the 
(Oc , Tc and Ac of the crowdion. In some sense, they can be considered as the actual 
parameters of the crowdion as they can be measured. Note that these parameters, 
as Ac are not well defined as they depend on the algorithm used to fit the numerical 
solutions. Figure [T2]represents a picture of c„ and a view of the phonon tail for 
similar to c„ where the perfect plane wave and its parameters can be appreciated. 


6 Excess energy and thermalized medium 

In this section we present the results of different simulations to show the capacity 
of the crowdions to survive a perturbed environment when larger energy is initially 
delivered and, second, the behaviour of the crowdions with temperature. 


6.1 Excess energy 

We present some examples of simulations when the lattice is given more energy 
than the 26.2 eV needed to produce the supersonic crowdion. The energies range 
from 130eV to 520 eV. They are represented in Fig.[T3j In (a) a single crowdion is 
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Fig. 13 (Color) Particle energy plots of several examples of crowdion formation in arbitrary units 
of ~ 3 eV. Initial energy increases form (a) to (f). Many features can be observed, among them the 
specific velocity of the crowdion V).. the formation of nonlinear waves and phonons, the formation 
of two crowdions and the survival of the crowdion in the severely perturbed media for hundreds of 
sites 
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formed after nonlinear waves are emitted. In (b) two crowdions are formed leaving 
behind an stationary linear wave. Note how the second crowdion survives to the tail 
of the first and the common velocity 14 of both. In (c) the excessive energy destroys 
the second crowdion which transforms into a highly localized nonlinear stationary 
wave. In (d) the second crowdion survives again, while in (e) it is again destroyed. 
Extensive phonon radiation and wandering kinks can be seen in the latter figure. In 
(f) a second crowdion survives for 150 sites in a highly perturbed media but it is 
finally pinned down. 


6.2 Thermalized medium 

An interesting question is whether the crowdion can travel trough a previously ther¬ 
malized medium. This is not only a question of general interest but particulary im¬ 
portant for mica muscovite. As it has been calculated in Ref. Tracks in mica: 50 
years later El, the recording process of tracks happens a few kilometers under¬ 
ground under large pressure and temperatures of 700-1000 K. Although much more 
work is necessary, the answer is positive. For comparison Fig. [14] (a-b) shows two 
simulations at 300 K and lOOOK in the system without substrate potential where 
the kink survives over hundreds of lattice sites. It is not really surprising as, if we 
compare the energy of the crowdion 26.2 eV with the mean thermal energy of a par¬ 
ticle ksT, the crowdion energy is 1000 and 300 times larger at 300 K and 1000 K, 
respectively. 

In the case of including the substrate potential, as shown in Fig. [T4| (c-d) for 
300 K and lOOOK respectively, the crowdion can also travel for hundreds of sites 
of the previously thermalized media. As it was studied in Ref. |0, the crowdion al¬ 
ways has finite kinetic energy, but the final total energy of the kink, Ej^, is always of 
the order of magnitude of the Peierls-Nabarro (PN) barrier. The equivalent kinetic 
energy equivalent for the thermalized media is 0.005 (0.013 eV) at 300K and 0.016 
(0.043eV) at lOOOK in normalized and physical units . These values are far be¬ 
low the energy difference between the PN barrier and the kink energy. However, in 
some simulations, for temperatures of lOOOK the thermalization is not completely 
achieved due to appearance of nonlinear waves instead of phonons. Therefore, lo¬ 
calized peaks of the background vibrations can interfere with the crowdion where, 
in some cases, it can be trapped leading to a highly localized nonlinear stationary 
perturbation. Figure[T4|(d) shows and example of this situation, where the crowdion 
is eventually trapped forming an interstitial defect. 

Thermal effects discussed in this section lead to different survival path lengths of 
the crowdions. If the hypothesis of crowdions propagating in mica muscovite is cor¬ 
rect, they might be related with some of the tracks observed in the mineral. Other 
feature of the presented simulations worth remarking on is that the high equiva¬ 
lent temperature of the nonlinear tail radiation of the crowdion is likely to favour a 
change of structure and the formation of tracks. 
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Fig. 14 (Color) Particle energy plots of two crowdions travelling in a previously thermalized 
medium at (a,c) 300 K and (b,d) 1000 K, top without and bottom with substrate potential. Color 
bars are in 101ogjQ(£) units 
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7 Recoil energy of 


If the hypothesis of quodons being vibrational entities of ions of potassium is cor¬ 
rect, the most likely source of energy is the recoil from because a) the energy 
will be given directly to the potassium ion K+ , b) the relative abundance and decay 
frequency of and c) because of the energies involved as explained below. 

The two most abundant isotopes of potassium are the stable and isotopes, 

with 93.7% and 6.7% abundance respectively. The next most abundant isotope is 
with a very long half life of 1.248 x 10^ years and abundance of 0.0117%. This 
isotope is the most important source of radioactivity for humans. 

As shown in Fig.[T5]and table|2] the nucleus experiences decay through dif¬ 
ferent branches with two daughter nuclei ^’^Ca and Ar 1412011 . The main parameters 
of the decay are /, the intensity of a given branch in % and Q, which is the difference 
between the rest masses of the parent and daughter atoms. The difference of mass 
between atoms is better tabulated than between nuclei. As the atoms are neutral the 
mass difference between nuclei has to take into account the difference in the number 
of electrons in the neutral atoms. The available energy will depend on the rest mass 
of the parent and daughter nuclei and other particles. It will be obtained below for 
each type of decay. 

The decay branches, j3 — and p -f involve the emission of an electron or a positron 
and a neutrino. The electron or positron velocities are such that they have to be 
treated relativistic ally, while the recoil velocity of the much heavier nuclei can be 


40k 

(1.2504(30) ■ 10® a) 



Fig. 15 Sketch showing the different decays and branching of Reproduced with permission 
from: Pradler, J., Singh, B., Yavin, I.: On an unverified nuclear decay and its role in the DAMA 
experiment. Phys. Lett. B 720(4-5), 399-.404 (2013). Copyright Creative Commons BY 3.0 
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described classically. We will suppose an electron to simplify the language, but a 
positron can be equally described in what follows. The maximum recoil energy of 
the nucleus is obtained when the neutrino gets no kinetic energy. The recoil energy 
is much smaller than the electron energy, so it can be neglected in the energy cal¬ 
culations while due to its large mass, it is essential for the momenta balance. The 
electron maximum energy is Ee = meC^+E, being E the available energy in the de¬ 
cay, and E^ = + PgC^, where pe is the momentum of the electron. Considering 

the parent nucleus at rest, the momentum of the nucleus is identical to the momen¬ 
tum of the electron piv = Pe = (^ /c)(£e —and the maximum nucleus recoil 
energy is given by E^ = plf/{2mi^). The decays always involve the emission of a 
neutrino and may include the emission of a photon, either 7 from the nucleus of X 
from the electron shell, although the latter have much smaller energy and momen¬ 
tum and will be of no importance for K+ recoil. The neutrino can be considered as 
a massless particle as its rest mass it known to be below 2.2eV/c^. Therefore for 
photons or neutrinos their energy is given by £v,y,x = pc. If only a photon or a neu¬ 
trino is emitted the recoil momentum p^ is equal to the momentum of the photon or 
neutrino and trivially Ef^ = P^/{2mf^). If there are only two daughter particles the 
recoil energy has a single value. 

Other data of interest are the ionization energies of K and of the daughter nuclei. 
If the recoil energy is larger than the ionization energy of the atoms that interact, 
it can be used to ionize an atom or ion and the energy cannot be transferred to the 
neighbours. The ionization energies of the daughter atoms from decay can be 
seen in table [T] A examination of the possible ionization processes is done in the 
following subsection. 


Table 1 Table of ionization energies (eV) of atoms involved in decay El 


Element 

I 

II 

III 

IV 

V 

Ar 

15.76 

27.63 

40.74 

59.81 

75.02 

K 

4.34 

31.63 

45.81 

60.91 

82.66 

Ca 

6.11 

11.87 

50.91 

67.27 

84.50 


As the lattice is formed by K+ , it is probable that the second ionization of K, 
31.6 eV, is an upper limit for crowdions or single row kink energies. 


7.1 decay branches 

Here, we analyze in detail the different decay branches. A summary if presented in 
table | 2 ] 

The decay branch that leads to ^’^Ca is: 





24 


JFR Archilla, YuA Kosevich, N Jimenez, VJ Sanchez-Morcillo and LM Garcia-Raffi 


j3 decay: 

Decay with emission of an electron. 

With/j3^ =89.25% and mass difference between atoms Q = 1311.07keV |l20|. 
As the Ca atom has an extra electron, discarding the electron binding energy of a 
few keV, the mass difference between nuclei is Q + and the energy available 

when emitting an electron is £ ~ 2 + meC^ — mgC^ ~ Q which will be shared be¬ 
tween the electron and the antineutrino emitted. Therefore, the maximum kinetic 
energy of the electron or endpoint is almost equal to Q. The daughter nuclei of 
^‘^Ca have a continuous distribution of energy with a maximum of = 42 eV at 
the endpoint corresponding to a velocity V — 14.4 km/s. 

The proton number increases by one, but the number of electrons does not 
change, therefore the daughter ion would be Ca++ with 50.6 eV third ioniza¬ 
tion energy. This is a likely origin of quodons for the decays with recoil energy 
smaller than the 31.6eV K second ionization energy. The recoils with larger en¬ 
ergy will be able to deliver up to 10.4eV after the first collision that could pro¬ 
duce breathers but not crowdions. 

The following processes have ^^'Ar as daughter nuclei being the difference between 
the atomic masses Q — 1504.69 keV. As the Ar atom has an electron less than K, 
discarding the electron binding energies the mass difference between nuclei is ~ 
Q — mgC^ and the energy available depends on the specific decay. 

ECl; Electron capture with decay to ^''Ar excited state and J radiation. 

With le = 10.55%, an electron from the shell is captured, therefore the available 
energy is E Q — mgC^ + nigC^ ~ Q. In this decay a monoenergetic neutrino of 
44keV is emitted with negligible recoil (26 meV) and the daughter nucleus is in 
an excited state. Thereafter, the excited nucleus decays to the ground state with 
the emission of a 1460keV 7 ray ||20|. The corresponding K+ recoil energy of 
‘^^Ar is Ejc ~ 29.2eV with velocity V = 12.0km/s. As this is a two body process 
Ei- has only slight variations due to interactions with the shell electrons. 

As no charge is emitted from the ion K+ , the daughter will also be a monovalent 
ion of Ar^, with 27.7 eV second ionization energy. So there is some probabil¬ 
ity that the first Ar^ collision with K+ will further ionize Ar^. The remaining 
energy 1.3 eV will not be enough to produce a kink but may produce a breather. 
ECIh-CE: Electron capture with decay to excited state and conversion electron. 
This is actually a subset of the previous decay, but with a probability 1=0.001 per 
100 decays, the 1460 keV yray emitted can interact with the shell and deliver the 
energy to an electron, called a conversion electron. Save for a few keV of binding 
energy the 7 energy is converted into kinetic energy of the electron, with a recoil 
for the ion of 49.7 eV and 15.6 km/s. This is the largest energy of all the recoils. 
As an electron has been emitted from the shell, the daughter ion will be Ar++ 
with 40.8 eV third ionization energy. This ionization and the 31.6 eV second one 
ofK+ are likely to occur. The remaining energies of 8.8 or 18 eV cannot produce 
a crowdion but will be able to produce breathers. 
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EC2; Electron capture with direct decay to ground state. 

With probability I = 0.2%, the energy available as in the decay above is £ ~ Q = 
1504.69. There is a direct decay to the ground state of ‘^^Ar after the capture of 
a shell electron and the emission of a monoenergetic neutrino that takes most of 
the energy available E ~ 1504.69keV minus the electron binding energy which 
is only a few keV EDEl. The recoil energy is 31.1 eV. The shell emits a 3 keV 
Auger electron when another electron of the shell occupies the vacancy left by 
the captured electron, however, this has a negligible recoil. 

The daughter nucleus has lost a positive unit charge but also the shell has lost two 
electrons, the captured one plus the Auger electron. Therefore the daughter ion 
will be Ar++, which has too little energy for further ionization of Ar++ or K+ 
which need 40.8 and 31.6 eV, respectively. Therefore, it would be a likely source 
of crowdions but difficult to distinguish from the j3 ^ recoil. 
decay with positron emission. 

With very low probability 7)3+ = 0.001%, the available energy is the mass dif¬ 
ference between nuclei minus the mass energy of the positron emitted, that is, 
E Q — irieC^ — meC^ = Q — 2meC^ = 483.7keV. The energy E is shared be¬ 
tween a neutrino, the emitted positron and the daughter nucleus. Therefore, 
the positrons have a continuum of energies with a maximum one or endpoint 
483.7keV ||4]|8][20l, which leads to the maximum recoil energy £4 ~ lOeV and 
velocity of 7 km/s. 

As the atomic number is decreased by one unit to Z — 1, the initial ion K+ has lost 
a positive unit charge, but there has been no change in the number of electrons, 
thus the daughter ion will be a neutral Ar interacting with short range forces with 
the neighbouring K+. The first ionization energy of Ar is 15.8 eV, so, actually, the 
Ar atom has less of the required energy for ionizing itself or for further ionization 
of K+ and will be able to keep the 10 eV energy. This seems too little to produce 


Table 2 Table of decays for 


Decay 

r 

ECl 

ECl+CE' 

EC22 

(3+ 

Intensity 

89.25% 

10.55% 

0.001% 

0.2% 

0.001% 

T (keV) 

1311.07 

1460 

1460 

1504.69 

483.7 

Emitted charged particle 

e“ 


e“ 

e“ 

e+ 

Recoil from 

v+e“ 

7 

e~ 

V 

V+e+ 

Max Recoil (eV) 

42 

29 2^ 

49.7“ 

31.1“ 

10 

Daugther ion (A=40) 

Ca++ 

Ar+ 

Ar++ 

Ar++ 

Ar 

Max V (Km/s) 

14.4 

12“ 

15.7“ 

12.2“ 

7 

Ionization of daughter (eV) 

50.6 

27.7 

40.8 

40.8 

15.8 

Aq (e) 

+ 1 

0 

+1 

+ 1 

-1 


* Subset of EC 1 when the gamma is delivered to a shell electron; Monocromatic 

^ Direct decay to Ar ground state, recoil from neutrino emission; 3 KeV Auger e“ 

EC: electron capture; CE: conversion electron; T: energy available excluding rest masses 
Ionization energy of K+ 31.6 eV 
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a kink but may produce a breather. Due to their positive charge, positrons leave 
tracks in mica muscovite 025II35I . 

A study with the correlation of positron tracks, thickness distribution of quodon 
tracks and other characteristics could make it possible to confirm the nature and 
characteristics of quodons. See Tracks in mica: 50years later for more details lIZTil . 


7,2 Secondary processes 


Electron-positron pair production: 

This is a secondary process after the 7 ray emission of 1460.82keV considered 
above jSl. It needs the interaction of the / ray with a nucleus, and the produced 
positron and electron can share the energy in any proportion. The maximum re¬ 
coil energy corresponds to a single particle taking almost all the energy except 
for the small amount taken by the nucleus, which is necessary due to momentum 
conservation. The available kinetic energy is £ = Ey — 2meC^ = 437.4keV and 
the maximum recoil energy is Ei^ = ^.^ eV. The probability of the combined pro¬ 
cess of electron capture and pair production is of the same order of magnitude as 
emission and also the energies are similar IS]. The probability of interaction 
of the / ray with a nucleus is proportional to which favors the interaction with 
potassium; however, potassium atoms are only 5% of the atoms in mica. 

As the energy is smaller than the second ionization energy of K of 31.6eV it is 
likely that the subsequent K+ -K+ collisions are elastic. 

Other secondary processes may also occur via other radioactive nuclei and their 
corresponding decay, but it will be beyond the objective of this work to continue the 
subject further. 


8 Summary 

In this work we have described in detail the magic mode for the strain or compres¬ 
sion of the bonds. A construction in terms of phasors has been developed in order 
to obtain an intuition of the relative phases and behaviour of the particles as the 
kink passes over them. We have considered an ID model for the close-packed lines 
of potassium ions inside a cation layer of mica muscovite using realistic potentials. 
There exists only a single kink with a specific velocity and energy dubbed the crow- 
dion. It is relatively well described by the magic mode but the kink is transformed 
into a double kink. It leaves behind a phonon wave with exponentially diminishing 
amplitude that travels at the same velocity as that of the kink. Simulations with dif¬ 
ferent initial energies bring about a variety of phenomena including the formation 
of two crowdions that leave behind nonlinear waves and phonons. The crowdions 
also survive at temperatures of 300-1000 K. Finally, an analysis of the possible de- 
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cay modes of has been performed including their possible consequences with 
respect to crowdion formation. A careful study of the tracks in mica muscovite com¬ 
pared with the decay modes can shed light on their characteristics and origin. 

The energy of the kinks or crowdions described in this work can be provided by 
the decay and is enough to expel an atom at the border. The crowdions survive 
to high temperature and travel long distances. They transport positive charge and 
therefore are very likely to be recorded in the form of dark tracks in mica muscovite. 
If they are the cause of the quodons or other marks observed in this mineral is still 
an open question. 
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